X-721 -66-77 




55 


NIKE TOMAHAWK 

PRELIMINARY PERFORMANCE STUDIES 


Q -w 0 
O JiJ 

"Mf 


cn 

N- 

:5k 


31 


N 


3 

vs 

Vs 

i 


GPO PRICE $ 
CFSTI PRICE(S) $ 


BY 

EDWARD E. MAYO 



FEBRUARY 1966 


«o» wuoi ajlihov* 



GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



X-721-66-77 


NIKE TOMAHAWK 

PRELIMINARY PERFORMANCE STUDIES 

by 

Edward E. Mayo 
Goddard Space Flight Center 
Greenbelt, Maryland 



NIKE TOMAHAWK 

PRELIMINARY PERFORMANCE STUDIES 


By 

Edward E . Mayo 
February 1966 


ABSTRACT 



In-house studies concerning the performance of the Nike Tomahawk 
are compiled. The following areas are investigated. 1) Nike Tomahawk 
particle trajectory characteristics. 2) Effect of payload weight on the 
Tomahawk mass, static margin, static stability, and natural frequency 
characteristics. 3) Tomahawk dynamic motions, and 4) Aerodynamic 
running load distribution over the Tomahawk body. The predictions 
should be updated as the vehicle properties and aerodynamic charac- 
teristics are refined. 
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UNITED STATES GOVERNMENT 

Memorandum 

TO : Flight Performance Section Files date: 15 March 1965 

FROM : Mr. Edward E. Mayo 

Flight Performance Section 

SUBJECT: NIKE- TOMAHAWK PARTICLE TRAJECTORY 

REFERENCE: (a) Stone* G. W. , and Connell* G* M. 2 Range Safety for the Nike-Tomahawk 

Rocket System With a Nine-Inch Diameter Payload* April 1964. 

(b) Letter of 21 January 1965 to Mr. John Lane from Sandia Corporation. 

Subject: Performance and Aerodynamic Data for the Nike-Tomahawk 

Rocket System. 

Nike-Tomahawk particle trajectory input values have been generated 
for the G.E. MASS Subprogram 1. The input data has been used to obtain 
trajectory data for a specific flight. The purpose of this memorandum 
is to document the Nike-Tomahawk thrust* vehicle weight and drag coeffi- 
cient histories for future reference. The values are given in proper 
increments for linear interpolation. Also* the trajectory characteristics 
for the specific flight are given. 

The Nike (M5E1) and Tomahawk (416) sea level thrust- time histories 
are given in Tables I and II* respectively* as obtained from reference 
(a). The Nike loaded weight (including interstage adapter) and pro- 
pellant weight were obtained from Mr, N. Peterson* Vehicles Section, as 
1320 and 750 lbs,* respectively. The Tomahawk (less gross payload) 
loaded weight and propellant weight were obtained from reference (b) as 
540 and 397 lbs., respectively. The weight-time histories were obtained 
by assuming that the weight of the vehicle is in proportion to the ratio 
of the impulse at a given time to the total impulse* that is* 


propellant 


The integrals were evaluated by the use of a planimeter. The resulting 
weight-time histories are given in Table III and Table IV. Vehicle drag 
coefficients for the nine-inch diameter payload are given in Tables V* 
VI* and VII as obtained from reference (a). 

1 
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Flight characteristics for a Wallops Island, 125 lb. nine-inch 
diameter payload, 15 foot length launcher, 80° elevation launch, are 
given in figure 1. The initial conditions for the flight are as 
follows: 


Run: 

Launch Site: 
Launcher Length: 
Longitude: 
Latitude: 

Launch Altitude: 
Launch Attitude: 
Launch Azimuth: 
Initial Velocity: 
Tomahawk Ignition 

The effects of varying gross 
on the Nike-Tomahawk system 
respectively. 


07 8 

Wallops Island 
15 ft. 

-75.480 deg. 

37.830 deg. 

12.81 ft. 

80.0 deg. 

-80.44 deg. (Inadvertent value) 

1.0 f ps 

16.0 sec. 

and Tomahawk ignition time 
given in figures 2 and 3, 


Time: 


payload weight 
performance are 


In summary, we now have the capability of simulating the Nike- 
Tomahawk flight with a particle trajectory,, Efforts are underway to 
determine the vehicle inertial and aerodynamic stability characteris 
tics as a function of payload weight. 


Edward E, Mayo 


Enclosures: 


(1) 

Table 

i. - 

(2) 

Table 

ii. - 

(3) 

Tc ble 

in. 

(4) 

Table 

IV. - 

(5) 

Table 

V. - 

(6) 

Table 

VI. - 

(7) 

Tab 1 e 

VII. 

(3) 

Figure 1. - 
Gross Paylo 
15 foot Lau 

(9) 

Figure 

2. - 


Nike (M5E1) Sea Level Thrust. 

Tomahawk (416) Sea Level Thrust. 

- Nike-Tomahawk Weight-Time History. 
Tomahawk Weight-Time History. 
Nike-Tomahawk Thrusting Drag Coefficient. 
Tomahawk Thrusting Drag Coefficient. 


80 deg. ; 


Tomahawk Apogee Altitude, Range at Impact, Time to Apogee 
and Time to Impact. Wallops Island; Launch Attitude = 80°; 
15 foot Launcher; Tomahawk Ignition Time = 16 sec. 
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(LO) Figure 3. - Effects of Tomahawk Ignition Time on the Nike- 
Tomahawk Apogee Altitude. Wallops Island; Gross Payload 
Weight = 125 lbs.; Launch Attitude = 80°; 15 foot Launcher. 

cc: Mr. K. R. Medrow 

EEM: skd 
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TABLE I 


Nike (M5E1) Sea Level Thrust* 


Time. Sec. [ 

1 

Thrust, lbs. 

1 

i 

0 

1860. 

.01 

25570. 

0 

O 

38121. 

.05 

40259. 

.09 

41003. 

.15 

41189. 

.84 

41840. 

1.14 

42119. 

1.74 

43235. 

2.04 

44165. 

2.34 : 

45280. 


1 

; Time. Sec, 

Thrust. lbs 0 

2.49 

45187. 

2.64 

44072. 

2.79 

41561. 

2.99 

36727. 

i 3.09 

32542. 

3.21 

23245. 

! 3.28 

! 15806. 

3.34 

| 10 507. 

3.40 

i 6323. 

3.46 

3161. 

3.54 

0 


* Nozzle exit area = 1.4630 sq. ft. 
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TABLE II 


Tomahawk (416) Sea Level Thrust.* 


Time* Sec* 

1 

Thrust, lbs. 

0 

0 

.10 

13847. 

.17 

12668, 

.22 

12766. 

.34 

11736. 

.90 

12030. 

1.10 

12080. 

1.60 

11785. 

2.20 

11588. 

3.00 

11736. 

3.50 

11686. 


Time, Sec, 

i 

Thrust, lbs. 

5.00 

11293. 

5. 30 

10999 8 

6.00 

9968 . 

7.00 

8888 o 

V ,H1 

8 347. 

8.00 

8004. 

8.1 *3 

80 53. 

8 c 30 

7268. 

8 0 7 2 

5892. 

9.20 

442. 

9 o 50 

0 


* Nozzle exit area = 0.4035 sq. ft. 
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TABLE III 


Nike-Tomahawk Weight-Time History.* 


Time, Sec. 

! Weight* lbs* 

0 

t 

1360.00 

.05 

1851.82 

.1 

1340.42 

.2 

1018.22 

1.2 

1589.10 

1.8 

1448.40 

2.4 

; 1302.15 

2.8 

1206.30 

3.0 

1165.00 

3. 15 

1 1140.00 

3.35 

1118.00 

3.45 

1112.00 

3.54 

1110.00 


* Less gross payload 


6 



TABLE IV, - Tomahawk Weight-Time History.* 


Time, Sec, j Weight, lbs, 

i 

540.00 


0 

.25 

.50 

4.00 

5.00 
6.00 
7.00 
8.00 
8.50 
8.70 
9.0 
9.5 

999. 


529.64 

517.29 

341.82 
292.79 
248.05 
207.99 
171.62 
155.307 
150.00 

143.83 

143.00 

143.00 


* Less gross payload 
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TABLE V. - Nike^Tomahawk Thrusting Drag Coefficient.* 


1 1 

Mach No, 

V 

' 0 

f 

; .360 

.50 

[ .365 

.57 

; .367 

.60 

j .370 

,64 

1 .380 

I .70 

: .400 

.73 

: .425 

1 .79 

: .47 5 

i .82 

.525 

.86 

.600 

.89 

.■702 

.96 

1 .951 

i .98 

: .990 


1 

Mach No. 

: s T ; 

; t 

1,00 

[ i 

.995 

1.06 

.990 1 

1.24 

.923 ; 

1.45 

.853 

1.70 

.77 5 . 

1.88 

.730 : 

2.06 

.690 , 

2.24 

.650 ; 

2.40 

.625 1 

2.53 

.615 

2.68 

i .605 

2.88 

; .590 

3.08 

.580 


* 1. Aerodynamic reference area = 1.474 sq. ft. 
2. Nine-inch diameter payload. 
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TABLE VI 


Tomahawk Thrusting Drag Coefficient.* 


Mach No* 

S T 

Mach No. 


1.25 

.955 

3.10 

.525 

1.95 

.800 

3.35 

.500 

2.00 

.725 

3.50 

.490 

2.20 

.675 

3.75 

.470 

2.45 

.625 

4,10 

.450 

2.70 

.585 

4.40 

.440 

2.75 

.570 

4.90 

.415 

2.90 

.550 

5,85 

.385 

3.00 

.540 

9.00 

.290 


* 1. Aerodynamic reference area = 0.4418 sq. ft. 
2. Nine-inch diameter payload. 
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TABLE VII. - Tomahawk Coasting Drag Coefficient.* 



■a 

1.25 

1.165 

1.60 

1,000 

1.90 

.900 

2.20 

.800 

2.50 

.730 

2.90 

.650 

! 3.25 

.600 

3.50 

.565 

3.75 

. 540 


f 

1 

! Mach No. 

C °c , 

i 

j A. 00 

.520 

4.30 

.500 

4.75 

.47 2 

5.25 

.445 

I 6.00 

.415 : 

; 7.00 

.380 j 

; 9.oo 

.315 j 

j 999. 

t 

; - 315 | 
1 


* 1. Aerodynamic reference area = 0*4418 sq. ft. 
2. Nine-inch diameter payload* 
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FIGURE 2. - The effects of gross payload weight 

on the Nike-Tomahawk apogee altitude, 
range at impact, time to apogee and 
t ime t o jjnpa c t , 
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EFFK-T OF PAYLOAD WEIGHT ON THE TOMAHAWK MASS, STATIC MARGIN, STATIC 
STABILITY, AND NATURAL FREQUENCY CHARACTERISTICS 


RE FEE j :\ T CE : 


(a) Memo of 3 June 1964, Mr. Mayo to Mr. Sorgnit, Subject: Weight, 

Center of Gravity, Pitch and Roll Moment of Inertia Determination 
Program . 

(b) Memo of 9 July 1964, Mr. Mayo to 671.2 Files, Subject: Static 

Stability and Natural Frequency Program. 

(c) Stone, C. W. , and Connell, G. M. : Range Safety for the Nike- 

Tomahawk Rocket System With a Nine-Inch Diameter Payload. April 
1964. 

(d) Letter of 21 January 1965 to Mr. John Lane from Sandia Corporation. 

Subject: Performance and Aerodynamic Data for the Nike-Tomahawk 

Rocket System. 

(e) Memo of 15 March 1965, Mr. Mayo to Flight Performance Section Files, 

Subject: Nike-Tomahawk Particle Trajectory. 


A brief study has been conducted to define the effect of payload 
weight on the Tomahawk mass, static margin, static stability, and 
natural frequency characteristics. The computer programs developed in 
references (a) and (b) were utilized in the study. The program input 
characteristics were determined through the use of references (c) 
through (e) . The study was conducted for the Aurora Nitehawk configu- 
ration shown in figures 1 and 2. A uniform density payload, as illus- 
trated by the cross hatched region in figure 2, was assumed. The purpose 
of this memorandum is to document the results of the study. 

Computations were performed for gross payload weights of 50, 75, 
100, 125, 150, 175 and 200 pounds. The vehicle aerodynamic characteris- 
tics were determined for launch angles of 75, 80 and 85°. The assump- 
tions used in determining the mass characteristics were: 

1. Propellant center-of-gravity location is constant. 

2. Nose weight distribution is conical as illustrated by 
the cross hatch in figure 2. 

3. Payload is of uniform density. 

4. Pitch and roll moments of inertia vary linearly with 
burning time. 
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Tomahawk weight and inertia input characteristics are presented in 
Table I. The resulting weight, center-of -gravity location, and pitch and 
roll moment of inertia characteristics are presented in figures 3 through 
6, respectively. Tomahawk aerodynamic input characteristics are presented 
in Table II. The resulting static margin, static stability, and natural 
frequency characteristics are presented in figures 10 through 12, respec- 
tively. The velocity, altitude, Mach number and dynamic pressure histories 
for the 100, 150 and 200 pound gross payload weights at launch angles of 
75, 80 and 85 degrees are presented in figures 7 through 9 as obtained from 
the GE MASS particle subprogram using the vehicle characteristics presented 
in reference (e) . The 50, 75, 125 and 175 pound gross payload weight tra- 
jectory characteristics are not presented in figures 7 through 9; however, 
they are on file within the Flight Performance Section for future reference. 

The characteristics are presented herein without discussion; however, 
several important points are: 

1. The vehicle characteristics are only as valid as the afore- 
mentioned assumptions. 

2. From stability considerations, a gross payload weight of 
around 80 pounds is the minimum flyable payload weight. 

This weight corresponds to a maximum peak altitude of 224 
statute miles or 360 kilometers. (Based on a nominal 80 
launch angle.) However, based on the recommended minimum 
static margin of 2 calibers, the maximum attainable alti- 
tude is about 200 statute miles. 

3. The variation of the static margin, static stability and 
natural frequency characteristics is insensitive to 
launch angle. 


Edward E . Mayo 


Enclosures: (4) 

(1) Symbols 

(2) List of Figures 

(3) Tables (2) 

(4) Figures (12) 

cc: Mr. K. R. Medrow 

Mr. G. E. MacVeigh 

Mr. E. E. Bissell 

Miss E. C. Pressly 
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SYMBOLS 


J N 0 


"Me 


xx 


yy 


t 

w 


x 

CO 


Normal force curve slope at ©< - 0, per radian 
Static stability parameter at<=* = 0, per radian 
Roll moment of inertia, slugs ft^ 

Pitch moment of inertia, slugs ft^ 

Time, sec. 

Weight, lbs. 

Longitudinal location measured from base, ft. 
Natural frequency, cps 


Subscripts 


b . o. 
c *g- 

c. p. 

i 


M.E.F. 


M.L.F. 
P . A. 
P.L. 


Burn out 

Cent er-of -gravity 
Cent er-of -pressure 
Ignition 

Motor empty -I- fins 
Motor loaded + fins 
Payload adapter junction 
Propellent lost 
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LIST OF FIGURES 


Figure 

1 Nike-Tomahawk sounding rocket system. All dimensions are 
in inches. 

2 Tomahawk TE-416 sounding rocket. All dimensions are in 
inches . 

3 Tomahawk weight-time history. 

4 Tomahawk center-of-gravity location. 

5 Tomahawk pitch and yaw moments of inertia. 

6 Tomahawk roll moment of inertia. 

7 Nike-Tomahawk particle trajectory. Wallops Island; 13 ft 
launcher; launch angle = 75°. 

(a) Gross payload weight = 100 lbs. 

(b) Gross payload weight = 150 lbs. 

(c) Gross payload weight = 200 lbs. 

8 Nike-Tomahawk particle trajectory. Wallops Island; 15 ft 
launcher; launch angle = 80°. 

(a) Gross payload weight = 100 lbs. 

(b) Gross payload weight = 150 lbs. 

(c) Gross payload weight = 200 lbs. 

9 Nike-Tomahawk particle trajectory. Wallops Island; 15 ft 
launcher; launch angle = 85°. 

(a) 

(b) 

(c) 

10 Tomahawk 

(a) 

(b) 

(c) 


Gross payload weight = 
Gross payload weight = 
Gross payload weight = 

static margin. Wallops 

Launch angle =75° 
Launch angle = 80° 
Launch angle = 85° 


100 lbs. 

150 lbs. 

200 lbs. 

Island; 15 ft. launcher. 
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LIST OF FIGURES (Continued) 


Figure 

11 


12 


Tomahawk static stability. Wallops Island; 15 ft. launcher. 

(a) Launch angle - 75° 

(b) Launch angle - 80° 

(c) Launch angle = 85° 

Tomahawk natural frequency. Wallops Island; 15 ft. launcher. 

(a) Launch angle * 75° 

(b) Launch angle = 80° 

(c) Launch angle = 85° 
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Table I. - Tomahawk Weight and Inertial Input Characteristics 


W M.E.F. 

- 

143 lbs. 

(lyy)^.E.F. 

= 63.2 slugs 

ft.' 

W M.L.F. 

S ! 

540 lbs. 

(lyy)M.L.F. 

= 207 slugs 

ft. 2 

^ Xc -g*) M< E.F. 

K 

3.46 ft. 

^xx^M.E.F. 

= 1.70 slugs 

ft. 

^ x c.g. M.L.F . 

ZZ 

5.88 ft. 

^xx^M.L.F. 

= 2.59 slugs 

ft. 

( x c.g.)p.L. 

= 

6.75 ft. 

ti 

= 16 sec. 


X P . A. 

= 

11.87 ft. 

t b . o . 

= 25.5 sec. 



t. sec. 

w p . L . » lbs - 

t. sec. 

Wp p , lb S • 

16.00 

0 

21.00 

247.2 

16.25 

10.36 

22.00 

292.0 

16.50 

22.71 

23.00 

332.0 

16.75 

35.69 

24.00 

368.4 

17.00 

48.36 

24.25 

376.3 

18.00 

98.66 

24.50 

384.7 

19.00 

148.5 

25.00 

396.2 

20.00 

198.2 

25.50 

397.0 
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*1. Aerodynamic reference length = 0.75 ft. 

2. Aerodynamic reference area = 0.4418 ft.^ 
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Fiqitl-t 7 - Continued. 
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(A) Gross payload vciqhf- loo Us. 


Figure, i - Ni ke-Tomak auik 'particle, 'trajectory . klallops Island 1 /-T -ft. launcher^ 
launch anqlc *■ to*. 
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F'3 Urt t - Con+inue-tf. 
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Gc) Gross payload loo tU. 


F '<^ure. &.- Concluded 
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Tmistt. 


(A) Gross payload weight - too lbs. 


Figure. 1- Alike- Toma hawk port: ole trajectory Wallops Island ; If ft launcher; 
launch angle = Tf* 
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Ih) Gross payload ivacyht* IfO Iks, 


Figure. Confinutd . 
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7W, set. 

(&■) l A unci, == 7T* 

Figure, 10.- Tomahavk sta+fa me rj Oa //ops Is /and. ) it 4f launder. 
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OPTIONAL. FONM NO- 10 
MAY 1111 KDITION 
OtA OKN. MO. NO. V 

UNITED STATES GOVERNMENT 

Memorandum 

: Flight Performance Section Files date: 9 February 1966 


from : Mr, Edward E. Mayo 

Flight Performance Section 

SUBJECT: TOMAHAWK DYNAMIC MOTION STUDY 


REFERENCE: 


(a) Nicolaides, Dr, John D.: Stability of Free Flight Missiles. October 

1964, 

(b) Memo of 15 March 1965, Mr, Edward E. Mayo to Flight Performance Section 

Files, Subject: Nike-Tomahawk Particle Trajectory 

(c) Stone, G. W., and Connell, G. M. : Range Safety for the Nike-Tomahawk 

Rocket System with a Nine-Inch Diameter Payload, April, 1964, 

(d) Memo of 29 June 1965, Mr, Edward E. Mayo to Flight Performance Section 

Files, Subject: Effect of Payload Weight on the Tomahawk Mass, Static 

Margin, Static Stability and Natural Frequency Characteristics 

(e) Memo of 10 September 1964, Mr, Edward E. Mayo to Flight Performance 

Section Files, Subject: Effects of Payload Weight and Length on the 

Cajun Mass, Static Margin, Static Stability, and Natural Frequency 
Characteristics. 

(f) Memo of 11 August 1964, Mr, Edward E. Mayo to Flight Performance Section 

Files, Subject: Effects of Payload Weight and Length on the Apache 

Mass, Static Margin, Static Stability, and Natural Frequency Character- 
istics, 

(g) Stone, George W,; The Magnus Instability of a Sounding Rocket, AIAA 
Paper No, 66-62, AIAA 3rd Aerospace Science Meeting, January 24-26, 
1966, 

(h) Price, D. A. , Jr., and Nelson, E, 0.: Final Report for Aerobee 350 

Lock-In Study, Lockheed Missiles and Space Company, 1965, 


INTRODUCTION 

The purpose of the six-degree-of -freedom simulation study disclosed 
herein is to investigate the Tomahawk vehicle dynamics at (1) booster 
separation, (2) resonance, and (3) post burnout conditions. The effects 
of fin misalignment, thrust misalignment and eccentricity, roll rate, 
static margin, level of pitch damping, and center-of- gravity offset are 
singularly investigated. The vehicle dynamics at booster separation 
and resonance are compared with equilibrium solutions of reference (a). 
Post burnout coning motions are currently under investigation in a sepa- 
rate study, 
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CONDITIONS OF STUDY 


The six-degree-of -freedom simulation time interval was from Nike booster 
burnout (t = 3*5 sec.) to approximately atmosphere exit time (t = 50 sec.). 
The initial conditions for position and velocity at 3.5 sec. were obtained 
from the 80-degree sea level launch particle trajectory runs of reference 
(b). The initial conditions are: 


| Payload Weight, 

lbs. 

i — ] 

. 100 

f 

125 

150 

I 

200 | 

i 

'Part. Traj. Run 

No. 

084 

078 

1 

085 

090 

i t t sec. 


3.5 

3.5 

3.5 

3.5 

Rsf * A G* de 8» 


-75.480 

-75.480 

-75.480 

-75.480 ] 

Ref. 0 G , deg. 


37.830 

37.830 

37.830 j 

37.830 1 

Ag* de 8» 


| -75.483 

-75.483 

-75.483 

-75.483 ; 

®G» de 8« 


37.83 

! 

37.83 

: 37.83 

37.83 j 

h, ft. 


4569.25 

| 4503.25 

! 4438.75 

4315.25 j 

I 

V E , fps 


2494.36 

2457.57 

2421.68 

2352.92 | 

r e* de §* 


78.589 

78.572 

1 

78.555 1 

t 

78.521 

Ae» de ® # 

i 

! 

j 

i 

90 

■ 1 

i 90 

; i 

1 1 

90 

i 

90 i 

— 1 


The initial body orientation and body rates were assumed to be 4-degree 
angle— of -at tack and 720 deg/sec. roll rate, respectively. 

A 0.2 degree angular thrust misalignment and 0.01 foot thrust eccen- 
tricity was assumed with the point of application at the base of the vehicle 
A 0.2 degree fin misalignment and a fin center-of-pressure location 7 inches* 
ahead of the base was assumed. The thrust and fin malalignments are consis- 
tent with those used in dispersion studies performed by Sandia Corporation 
(reference c). In the study presented herein, the malalignments are always 
in the worst possible orientation. The vehicle mass and aerodynamic charac- 

e j 1 S j j? S . are presented in Tables 1 through XV as obtained via reference (d) 
and additional calculations. 

The six-degree-of-freedom simulation cases are summarized in the follow- 
ing table. 
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Run 

Payload j 

Weight, lbs. 

Thrust 

Malalignment 

Fin 

Misalignment 

Fin Cant* 
Min* 

m 

Kj* 

192 1 

100 

Yes 

Yes 

16 

16 

1 

187 

125 

Yes 

Yes 

16 

16 

i 

193 

200 

Yes 

Yes 

16 

16 

i 

190 

125 

Yes 

Yes 

18 

16 

i j 

I i 

191 

125 

Yes 

Yes 

14 

16 

i ! 

193 

125 

Yes 

Yes 

16 

16 

0 ! 

188 

125 

No 

Yes 

: 16 

16 

1 j 

189 

125 

Yes 

No 

16 

16 

1 

1 

194 

I 125 

i 

Yes 

Yes 

i 

; 16 
j 

12 

l 

4 

i 1 

j 1 


* G m q multiplier 


DYNAMICS 

Roll-Fitch RateS e - The Tomahawk vehicle has four (4) wedge-slab type fins 
which may be canted from 14 ! to 18' to produce the desirable roll history. 

The estimated roll forcing and roll damping coefficients are given in tables 
13 and 14* respectively* For canted fins* the roll rate is proportional to 
the vehicle velocity* An examination of the Tomahawk velocity-time history 
will reveal that canted fins are desirable as roll producing devices* This 
stems from (1) allowing the roll rate to remain safely below the pitch fre- 
quency during coast* (2) a rapid increase in roll rate through resonance 
during which the roll forcing moment is much greater than the damping moment* 
and (3) the large velocity increase during burning forces the roll rate to a 
safe distance above the pitch frequency* 

The effects of gross payload weight, Tomahawk ignition time and fin cant 
on the pitch-roll histories are given in figures 1 and 2. From figures 1 and 
2* the undamped pitch frequency prior to resonance is insensitive to payload 
weight* This was also observed in references (e) and (f) for the Capache 
vehicle. The primary effect of payload weight is the change in space roll 
rate which results from the change in burnout velocity* For the range of 
payload weights investigated (100 to 200 lbs.)* the space roll rates are 
within 5*5 + 0*5 cps* From figures 1(b)(1) and 1(b)(2)* there is essentially 
no effect of Tomahawk ignition time on final roll rate* 
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The effect of fin cant on the roll rate history is given in figure 2 
for a nominal payload weight of 125 lbs. and a Tomahawk ignition time of 
16 seconds. Also shown on figure 2 are equilibrium roll rates computed by 



( 1 ) 


Initially^ et t 3o5 sec 0 , the actual roll rate value is the assumed 2 cps 
initial condition. The roll rate reaches the steady state value within 
about 0.2 sec. and remains nearly equal to the steady state value until 
Tomahawk ignition. There is considerable lag during thrusting which is 
attributed to vehicle longitudinal acceleration coupled with roll inertia 
effects. After burnout, the actual roll rate overshoots the equilibrium 
value by approximately 0.5 cps at t = 50 sec. The space roll rate is pro- 
portional to the fin cant and the roll rate for all fin cants investigated 
is within 5.6 + 0.75 cps. 


Booster — Separ ation . - The initial disturbance at booster separation damps 
from 4 to 1° in approximately 1.1 second. See figures 3 and 4. This is 
in agreement with the 0.57 second Nutation and Precession Arm half life 
time computed from the equilibrium solutions given in reference (a). Hence, 
the disturbances of the Tomahawk at booster separation are heavily damped 
and should not impose any serious flight abnormality. 


R g s onance . - The effects of gross payload weight and fin cant on the com- 
bined angle-of-attack build-up at resonance are given in figures 3 and 4 
respectively; and the maximum combined angle-of-attack is compared with * 
equilibrium solutions of reference (a) in figures 5 and 6. The maximum 
combined angle-of-attack is approximately 4 degrees and is within 1 degree 
of the equilibrium value. From figure 5, the effect of varying payload 
weight (and, thus, minimum static margin) has a relatively small effect on 
the maximum combined angle-of-attack. This can be explained by referring 
to the static margin histories given in reference (d). The resonance time 
for the Tomahawk is approximately 20 seconds; whereas, minimum static mar- 

fon ?K CUr e C u 5 “ 5 SeC * ) * In goin S from a B ross payload weight of 

200 lbs. to 100 lbs., the minimum static margin decreases from 4 calibers 

to 1 caliber; whereas, at resonance time, the static margin decreases from 
5.25 to approximately 3.25 calibers. Similarly, for the case of zero mini- 

static M margin of over 2 calibers exists at the resonance 
time. The small change in f[ max with varying fin cant exhibited in figure 6 
is expected since, referring to figure 2, the resonance time spread for ex- 
treme fin cants is only 1.25 seconds. 
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The effects of the various trim producing malalignments and aerodynamic 
parameters on the vehicle dynamics were singularly investigated. The effects 
of thrust malalignment, fin misalignment and level of pitch damping are given 
in figure 7 and the vehicular constraints, together with level of pitch damp- 
ing and axial force contributions to the equilibrium resonance build-up are 
presented in the following table* 




' 

b 

= 16', W p 

= 125 lbs., t 2 

= 16 

sec* 




! 

| ! 

I 

p ! 

CJ 

1 

°2 



m 

■ 

K-, 

J Q 

K 

K 3 

1 

K 3 

Motion 

: l 

rad/ sec 

rad/ sec 

rad/ sec 

1/ sec 

. 

1/ sec 

■ sec 

i 

sec ■ 

... _ _ j 

deg 

K 3 

J 0 

deg 

] 

P & Y ! 

14.89 

14.43 

-14.33 

-.37 

-.37 

•1.86 

j 

1.87 

.60 

18.95 

11.37 ! 

P, Y, H 
& S* 

14.89 

14.43 

-14.33 

-.60 

-.60 

| 

1 1. 15 

1.15 

.60 

11.75 

7.05 . 

, P» Y, H, ! 
S W & S u 

14.89 

14*43 

-14.33 

-.82 

-.82 

! .85 

| 

.85 

.60 

8.67 

i 

5*20 

; y, h, 

Swt ^u» 
St C m ^=0 

14.89 

j 

j 

j 14.43 

-14.33 

-.53 

i 

j 

-.53 

1 

1 

i 

il.32 

1.32 

j 

.60 

! ! 

i 

13.47 

1 

i 

I | 

| 8*08 | 

I | 

: Y, H, 

Sw* S us 
& C x = 0 

i 

| 

| 14.89 

| 

i 

l 

! 14.43 

i 

j 

' 

-14.33 

-.83 

1-.83 

i 

j .83 
j 

| .83 

i 

*60 

8.49 

5.09 

Where 

P & Y 


Pitching 

• 

and Yawing 




1 


i 

P, Y, H, 

^W 

Pitching 

, Yawing, Heaving and Swervin 

§ 


f 

P, Y, H, 

S ws s u 

Pitching. 

, Yawing, Heaving, 

Swerv] 

ng and Surging 


i 


From figures 7(a) and 7(b), the fin misalignment accounts for approximately 
40 percent of ^ max , while the thrust malalignment accounts for the remaining 
60 percent. The effect of level of pitch damping on the motion is given in 
figure 7(c). From figure 7(c) and 6, a 15 percent increase in ^ max occurs 
if ass umed zero* This is considerably less than the 55 percent increase 
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predicted from equilibrium solutions. See above table. The level of 
pitch damping had a pronounced effect on the post burnout coning motion. 

As seen from figure 8, the Nutation and Precession Arm damping rates 
after burnout decrease exponentially with time, (Figure 8 is for full 
estimated C m , ) Hence, the degradation in damping due to decreasing 
density, aerddynamic coefficient degradation and post resonance residual 
motions can result in large post burnout coning motions. It should be 
noted that the post burnout coning build-up exhibited in I inure 7(c) 
serves only as an indicator that the damping rate of the Nutation and/or 
Precession Arm is sour. In the study herein, the Magnus term (which ap- 
pears with Cmq in determining the Nutation and Precession Arm damping 
rates) was neglected. In reference (g), the Magnus moment coefficients 
were extracted from divergent flight characteristics of the Tomahawk vehi- 
cle and were shown to be large, highly non-linear and to vary rapidly with 
Reynolds number. The post burnout coning motion of sounding rockets is 
currently under investigation in a separate study by Mr, James McGarvey, 
a Fairchi ld-Hil ler employee. 

Lateral Center-of -Gravity Limit . - As noted previously, the Nutation and 
Precession Arms are heavily damped prior to resonance, and y\ max at reso- 
nance is in agreement with equilibrium values predicted by resonance 
instability theory. Thus, ^ raax can be estimated by procedures such as 
described in the Appendix. The resulting max versus time is given in 
figure 9 (p = O ). Figure 9 illustrates that if "lock-in" occurs, max 
will grow exponentially until flight failure. The allowable center-of- 
gravity offset (or induced rolling moment) to insure a break-out of 
"lock-in" was estimated according to the procedures of reference (h). 

Considering 4c.G, = 0, then the induced rolling moment required 
to maintain "lock-in 11 may be approximated by 

- Cl! - S C lj + C lp (<i|) (2) 

The resulting Ci values are given in figure 10 as a function of time 
and fin cant angle. Figure 10 is interpreted as follows. If an induced 
rolling moment of -0.015 exists, then for & = 16' roll "lock-in" will 
be maintained until 20 sec. ( ^ max 6° from figure 9) at which time 
break-out occurs. If for the same cant angle, an induced rolling moment 
coefficient of -0.030 exists, then roll "lock-in" will be maintained 
until 30.5 seconds ( max = 20° from figure 9) at which time break-out 
occurs. From equation 2 and as seen in figure 10, the allowable induced 
moment is a direct function of the fin cant angle. 
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Considering C]^ = 0 and the worst center-of-gravity orientation, the 
allowable center-of-gravity offset may be approximated by 


4 C.G. 


12d 


( S c. 


+ c 


max 


1 P ) 


(3) 


The resulting £C,G, values are given in figure 11 as a function of time and 
fin cant angle. Figure 11 is interpreted similar to figure 10, For f ■ 16', 
a center-of-gravity offset greater than 0,165 inches can prevent break-out 
from "lock-in"; whereas, for 4 C.G. X0.11 inches, break-out occurs at 
t = 20 sec. ( ^ max = 6° from figure 9). From equation (3) and as seen in 
figure 11, the allowable center-of-gravity offset is a direct function of 
the fin cant angle. 


CONCLUDING REMARKS 

The significant conclusions of the enclosed study may be summarized 
as follows: 

1* From consideration of the vehicle's inherent velocity- time his- 
tory, the canted fin roll control design appears attractive. 

2, The undamped pitch frequency prior to resonance is insensitive 
to payload weight. 

3* The space roll rates for all payload weights and fin cants in- 
vestigated are within 5*6 + 0.75 cps. 

4* The disturbances at booster separation are heavily damped and 
should not impose any serious flight abnormality. 

5* The maximum combined angle-of-at tack at resonance is approximately 
4 degrees and is within 1° of the equilibrium value. The build-up 
at resonance is insensitive to payload weight and fin cant. 

6. Fin misalignments account for approximately 40 percent of the com- 
bined angle-of-attack build-up at resonance while thrust malalign- 
ments account for the remaining 60 percent. (Only fin and thrust 
malalignments were considered.) The build-up is insensitive to 
the level of pitch damping. 
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7# The degradation in damping due to decreasing density, aerodynamic 
coefficient degradation and post resonance residual motions can 
result in large post burnout coning motions. 

8, Considering i * 16', 4 C f G, * 0, an induced rolling moment 

coefficient, Cj^ . # of -0,030 is sufficient to sustain roll "lock- 
in" until t = 30.5 seconds ( >\ = 20°). Considering £ = 16', 

c li = 0, a center-of-gravity offset, A C.G., greater than 0.165 
inches can prevent break-out from roll "lock-in" and, thus, flight 
failure. 
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APPENDIX 


RESONANCE INSTABILITY THEORY 


As evidenced by the six-degree-of-f reedom runs and substantiated from 
equilibrium solutions, the Nutation and Precession Arms arc heavily damped 
( 0*6 sec*) and reach extremely small magnitudes prior to resonance. 

Hence, the purpose of the equilibrium solutions herein was to establish 
how well the combined angle-of-attack at resonance could be predicted by 
resonance instability theory* 

From reference (a), the rolling trim arm, K3, at resonance is given 
by 


K 



*2 - i *1 Cp~- v 2 ) 


(Al) 


and the non-rolling trim arm, 1(3, is given by 


k 3 


e 


n 3 

N^ 


(A2) 


Dividing (Al) by (A2) and rationalizing, the magnification factor K3/K3 
may be expressed as 


Ko/K- 


Ai(A 2 ) 2 + A l Cp - ^ 2 ) 




(p 


- 0 2 )2 


(A 3 ) 


Neglecting Magnus effects, 


No = - a*. 

2 1 

Where = C moc qsd 

Neglecting Magnus and apparent mass, the damping rates, and^^ 2 , for 

the Nutation and Precession Arms are given by 
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c 2 qs 


M r 


2Vm ti--e, + a (ltT , 




- Cg °^ qS (i +r ) + ci .r ) 

2vm 21 


Where M n = M n + 

q q Aero <ljet 


" C m q ^>qsd + m(x cg - x e ) 2 - I 



The rotation rates, u l and ^2* f° r t ^ ie Nutation and Precession Arms are 
given by 


£is- (i - 

21 ^ 


0 2 


P I 


21 


X (1 + 1 ) 


The non-rolling trim arm K3 is given by 


T(x cg - x e )(sin * 


M ( 


21 


m 
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and the magnitude of the rolling trim at resonance m^y be determined from 
The Nutation and Precession half lives are given by 



Y~ In ( 1 / 2) 

A 2 


The foregoing equations were programmed by C. Hutton to yield the 
equilibrium solutions enclosed herein. 


61 


SYMBOLS 


C D 

\ 

C N* 

C m, 

Cjn 


Drag coefficient 


ij J 

Induced rolling moment coefficient, — — 

qsd 

Normal force coefficient curve slope at 


Pitching moment curve slope at = 0 

^ ^m 


Damping coefficient. 





- c D 

' C NoC 

Aerodynamic reference length, *75 ft 0 
Nutation arm half life 

Precession arm half life 
Altitude 

Pitch moment of inertia 
Roll moment of inertia 

C multiplier 

q 

Rolling trim 
Non-rolling trim 

Vehicle mass 
Induced rolling moment 
Mach number 
Roll rate 
Pitch rate 
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SYMBOLS 


£ 

T 

t 

fc 2 

V 

w 

x 

7 

oe 

r 

4 C.G. 

$ 

if 

£ 

A 2 

A 

r 


o 

Aerodynamic reference area, 0.4418 ft* * 

Gyroscopic stability factor 

Thrust 

Time 

Second stage ignition time, sec* 

Velocity 

Weight 

Longitudinal distance measured forward of base 

Lateral distance measured from x-axis 

Angle-of -attack 

Flight path angle 

Center-of-gravity lateral offset 

Fin cant angle 

Fin misalignment angle in X-Z plane 
Thrust misalignment angle 
Combined angle-of -attack 
Nutation damping rate 

Precession damping rate 
Geodetic longitude 
Azimuth 



s 
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SYMBOLS 


U 


1 


Nutation arm rotation rate 


OJ 


2 


Precession arm rotation rate 


SUBSCRIPTS : 

C Coasting 

c # g. Center- of -gravity 

Cop* Center-of -pressure 

e Nozzle exit 

E Relative to the Earth 

f Fin 

i Induced 

max Maximum 

p Payload 

t PDint of application of thrust 

T Thrusting 
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Figure 


I 


2 

3 


4 


Effects of gross payload weight on the Tomahawk undamped pitch 
rate and roll rate histories. 


(a) Gross payload weight = 100 lbs,, t 2 = 16 sec, 

(b) Gross payload weight = 125 lbs, 

(1) t 2 = 16 sec. 


(2) t 2 = 12 sec, 

(c) Gross payload weight = 200 lbs,, t 2 = 16 sec. 

Effect of fin cant on the Tomahawk roll rate history. Gross 
payload weight = 125 lbs,, t 2 = 16 sec. 


Effect of gross payload weight on the Tomahawk combined angle- 
of-attack history, $ - 16', £ = 0,2°, J . = ,01 ft,, 

if = -0.2 deg. 0 ^ 


(a) 

Gross payload weight = 100 lbs*, t 2 = 16 sec. 

(b) 

Gross payload weight = 125 lbs. 
(I) t 2 = 16 sec. 


(2) t 2 = 12 sec. 

(c) 

Gross payload weight = 200 lbs., t 2 = 16 sec. 

Effect of 
history. 

fin cant on the Tomahawk combined angle-of -at tack 
£ = 0.2 deg., p t = .01 ft., $ f = -0.2 deg. 

(a) 

S = 14' 

(b) 

$= 16' 

(c) 

S = 18' 
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Figure 


5 


6 

7 


8 


Effect of gross payload weight on the Tomahawk maximum combined 
angle-of-attack at resonance, 6 = 16 ', 6 = 0,2°, 1 t = .01, 

S f = “0*2 deg. 

Effect of fin cant on the Tomahawk maximum combined angle-of- 
attack at resonance, £ = 0.2°, J t = .01 ft., $ £ = "0.2 deg. 

Effect of thrust malalignment, fin misalignment, and level of 
pitch damping on the Tomahawk combined angle-of-attack history. 

& = 16', e = 0.2°, Jt = .01 ft., S f = -0.2°, Kj_ = 1.0, 
t 2 = 16 sec. except as noted. 

(a) 6-0 deg., = 0 ft. 

(b) $ f = 0 deg. 

(c) = 0 C m ^ = 0 

Equilibrium Nutation and Precession Arm half life history. 


9 Equilibrium combined angle-of-attack history for M locked-in M 

condition. 6 = 0 e 2 deg,, p = ©01 ft., $£ = -0.2°. 

10 Magnitude of induced rolling moment coefficient required to sus- 
tain "locked-in" condition. £ = 0.2°, 2t~ 

Si = -0.2°, A C.G. = 0. * 


11 Magnitude of center- of -gravity offset required to sustain 

"locked-in" condition, £. = 0.2°, 1. = .01 ft., £ f = -0.2°, 

G li = o. d 


67 



TABLE I 


TOMAHAWK COASTING DRAG COEFFICIENT* 


M 

1.25 

1.6 

1.9 

2.2 

2.5 

2.9 

3.25 

3.5 

3.75 
4. 
4.3 

4.75 

5.25 

6 . 

7. 

9. 

999. 


1.165 

1 . 

.9 

.8 

.73 

.65 

.6 

.565 

.54 

.52 

.5 

.472 

.445 

.415 

.38 

.315 

.315 


^Aerodynamic reference area 


0.4418 ft 2 
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TABLE II 


TOMAHAWK THRUSTING DRAG COEFFICIENT* 


M 

c d t 

1.25 

.955 

1.95 

.8 

2. 

.725 

2.2 

.675 

2.45 

.625 

2.7 

.585 

2.75 

.57 

2.9 

.55 

3. 

.54 

3.1 

.525 

3.35 

.5 

3.5 

.49 

3.75 

.47 

4.1 

.45 

4.4 

.44 

4.9 

.415 

5.85 

.385 

9.0 

.29 


O 

★Aerodynamic reference area = 0.4418 ft^ 
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TABLE 111 


TOMAHAWK NORMAL FORCE COEFFICIENT CURVE SLOPE* 


M 

Cm 

H oc 

m ; 

CN ~ 

1*2 

38.96 

4.8 

10.43 

1.43 

23.19 

5.15 

10.08 

1.6 

23.66 

5.6 

9.74 

1.8 

22.0 

6.0 

9.51 

:\o 

19.6 

6.55 

9.28 

2.23 

17.88 

7.0 

9.11 

2.5 

16.44 

7.5 

9.0 

3.0 

14.1 

8.0 ! 

I 

8.94 

3.4 

12.95 

1 

8.5 

8.82 

3.75 

12.03 

8.85 

8.82 

4.1 

11.34 

10.0 

8.82 

4.45 

10.77 




* 1 * per radian 

2. Aerodynamic reference area = 0.4418 ft^ 
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TOMAHAWK PITCH DAMPING COEFFICIENT 
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** Current time 


TOMAHAWK PITCHING MOMENT COEFFICIENT CURVE SLOPE 
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per radian 2 # Aerodynamic reference area = 0.4418 ft 2 
Aerodynamic reference length - 0*75 ft 



TABLE VI 


TOMAHAWK SEA-LEVEL THRUST 


t. sec* 

Thrust, lbs. 

t 

t. sec* 

Thrust, lbs. 

0 i 

0 

5. 

11293 

.1 j 

13847 

5.3 

10999 

.17 

12668 

6. 

9968 

.22 

12766 

7. 

8888 

.34 

11736 

7.81 

8347 

.9 

12030 

8. 

8004 

1.1 

12080 

8.15 

8053 

1.6 

11785 

8.5 

7268 

2.2 

11588 

8.72 

5892 

3. 

11736 

9.2 

442 

3.5 

11686 

9.5 

0 


* Phase time 
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TOMAHAWK CENTER -OF -GRAVITY LOCATION 



74 


TOMAHAWK GROSS WEIGHT, W, LBS 



75 


Current time 



TOMAHAWK ROLL MOMENT OF INERTIA, I x , SLUGS FT 
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TOMAHAWK PITCH AND YAW MOMENTS OF INERTIA, I, SLUGS FT 2 




TABLE XI 


TOMAHAWK CENTER-OF -PRESSURE LOCATION* 


1 

M 

x cp* 

1.2 

1.523 

2.3 

3.192 

2.75 

3.718 

3.0 

4.002 

3.5 

4. 468 

4.0 

4.865 

4.5 

5.218 

5.0 

5.518 

5.8 

5.93 

6.4 

6.23 

7.0 

6.51 

8.0 

6.905 

9.0 

7.258 

10.0 

7.258 


* Measured forward from base 




TABLE XII 


TOMAHAWK FIN LIFT COEFFICIENT* 


M 

C k®c 

F 

M 

Cl 

L OC 

F 

0 

15.563 

2.5 

10.815 

.2 

15.659 

2.75 

9.748 

i 

.4 

15.958 

3.0 

8.924 

• 6 

16.51 

3.5 

7.583 

.8 

17.42 

4.0 

6.567 

1.0 

18.946 

4.5 

5.842 

1.2 

20.555 

5.0 

5.252 

1.4 

18.899 

6.0 

4.382 

1.6 

17.266 

7.0 

3.742 

1.8 

15.369 

8.0 

3.266 

2.0 

13.814 

10.0 

2.606 

2.25 

12.144 

15.0 

1.732 


* 1. per radian 

2 

2* Aerodynamic reference area = 0.4418 ft 
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TABLE XIII 


TOMAHAWK ROLL FORCING COEFFICIENT* 


M 

Cl < 

i 

M 

Cl # 

0 

31.977 

2.5 

22.221 

.2 

32.174 

2.75 

20.028 

.4 

32.789 

3.0 

18.335 

.6 

33.922 

3.5 

15.58 

.8 

35.792 

4.0 

13.492 

1.0 

38.926 

4.5 

12.004 

1.2 

42.233 

5.0 

10.79 

1.4 

38.831 

6.0 

9.003 

1.6 

35.476 

7.0 

7.689 

1.8 

31.578 

8.0 

6.711 

2.0 

28.382 

10.0 

5.354 

2.25 

24.952 

15.0 

3.558 


* 1. per radian 

2 

2. Aerodynamic reference area = 0.4418 

3. Aerodynamic reference length = 0.75 ft. 
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TABLE XIV 


TOMAHAWK ROLL DAMPING COEFFICIENT* 


* 


M 

..... _ . . 

Ci 

V 

M 

V 

0 

-34,805 

2.5 

-58,931 

.2 

-85.327 

2.75 

-53.11ft 

.4 

-86.959 

3.0 

-48.626 

.6 

-89.963 

3.5 

-41.318 

.8 

-94.924 

4.0 

-35.782 

1.0 

-103.236 

4.5 

-31.835 

1.2 

-112.004 

5.0 

-28.617 

1.4 

-102.982 

6.0 

-23.877 

1*6 

- 94.085 

O 

O 

r- 

-20.391 

1.8 

- 83.748 

8.0 

-17.798 

2.0 

-75.272 

10.0 

-14.198 

2.25 

- 66.175 

15.0 

- 9.437 

i 


■) c 4 


2* per radian 

2 

3* Aerodynamic reference area = 0*4418 ft 
4 # Aerodynamic reference length = *75 ft 
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TOMAHAWK JET DAMPING 





5 . 10 . 15 . 20 . 25 . 90 . 96 * 40 . 45 . 50 . 

TIME FROM LRUNCH. SEC. 
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NIKE TOMRHRUK RUN (81 



10 . 


15 . 


20 . 25 . 80 . 

TIME FROM LAUNCH. SEC. 


86 * 


to. 


t5. 


50 . 
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NIKE T0MRHRWK RUN 193 



5 . 10 . 15 . ZO* Zf. 30 3 sT HC 


TIME FROM LRUNCH. SECo 
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NIKE TOMHHHWK RUN 192 
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TIME FROM LAUNCH. SEC. 
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NIKE TOMRHRRK RUN 193 
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6. 10. 16. 20. 26. SO. SE. *0. *£. 50. 

TIME FROM LAUNCH. SEC. 
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20 . 26 . 30 . 

TIME FROM LAUNCH. SEC. 
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COMB I 



5. 10. IS. 20. 25. 80. 85. 40. 15. 50. 

TIME FROM LRUNCH, SEC. 
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UNITED STATES GOVERNMENT 

Memorandum, 

TO : Flight Performance Section Files date: 25 August 1965 

FROM : Mr. Edward E, Mayo 

Flight Performance Section 

SUBJECT: APPROXIMATE SECOND ORDER SHOCK EXPANSION METHOD 


REFERENCE: (a) 
(b) 


(c) 


(d) 


Memorandum of 23 March 1965, Mr. E. E. Mayo to Flight Performance 
Section Files, Subject: Capache Running Load Distribution 

Syverston, Clarence A. and Dennis, David H.: A Second-Order 

Shock-Expansion Method Applicable to Bodies of Revolution Near 
Zero Lift. NACA Report 1328, 1957. 

Memorandum of 7 July 1965, Mr. E. E. Mayo to Flight Performance 
Section Files, Subject: Tangent Ogive Geometric and MASS Char- 

acteristic Equations 

Ames Research Staff: Equations, Tables and Charts for Compressible 

Flow. NASA Report 1135 


INTRODUCTION 


In reference (a), a method is given whereby the running load distri- 
bution for a cone-cylinder configuration may be calculated. For ogive- 
cylinder configurations, the second-order shock expansion method of 
reference (b) is the best practicle theory for determining the running 
load distribution. At present, a computer program is being developed 
for the second— order shock expansion method in which the program input 
will consist of the Mach number, body coordinates and tangency points. 

To satisfy current needs until the program is developed, the pro- 
cedures for application of the approximate second-order expansion theory 
have been established. The purpose of this memorandum is to outline 
the procedures for determining the running load distribution for both 
the cone-cylinder and ogive-cylinder configuration via the approximate 
second-order shock expansion method given in reference (b) . As an ex- 
ample of application of the procedures, the running load for a cone- 
cylinder and ogive-cylinder configuration are presented in the appendixes. 
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ANALYSIS 


Cone-cylinder configuration . - The nondimen sional running load distribution 
for the cone portion may be determined by equation ( 5 ) of reference (a) # ie # 


_ n 

ocqd 


c 


tanr Cn <* c 


( 1 ) 


The running load distribution over the cylindrical portion may be via 
equation (Cl) of reference (b). From equation (Cl) 

-G2 fn -G2fa 

C N | = G 1 e (1 “ e ) 

°c la 


Differentiating with respect to x/d yields 

d C]\j ^ I r __ ~| 

— 3 / = G,G o exp / -G (fn + x/d) I 

d(x/d) j a 2 L 1 J 

The relationship between n/*c qd | a and d C Nocl 

d(x/d) | 

rr dC N«^ 


is given by 


n/oc qdj a - — 


4 d(x/d) 


Hence, the running load distribution over the cylindrical afterbody is 
given by 


n 

<*qd 


a 


rr 

4 


G 1 G 2 



G 0 (fn + x/d) 


il 


( 2 ) 


The use of equations ( 1 ) and (2) will yield the running load distribution 
over the cone-cylinder configuration. 


Ogive-cylinder configuration . - The running load distribution for the ogive 
is given by equations (C6), (C 7 ) and (C8) of reference (b) as 


-A- = 


j tan £ C»j I + G, tan £ + £/l - Gq) C Noc/ tcv - tan 2 $ 

3 * V j tcv 4 ° tan & v 


( 3 ) 
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where 


G _ £fa 8 ln 2 A, 

3 ~ P v sin 2.4 

G 4 = 2(1 - e ^ a ) 


- ^ 


K is given by eq # (C5) of reference (b) as 

i, = t ^ Ps a /Po ° ) Ms a ifkv 

r “ 2(1 - P s 7P«,)(M 2 - 1) Ji 


Where -j is the one-dimensional area ratio given by 

v +i 

i )Lzl 2 I 2( V-1) 

; _ i I 1+2 


0 M y +1 


From eq* 14 $ reference (b) 


CNoc = r d x 


Which, upon nondimensionalizing, yields 

r 1/d 

C’Njp = 8 \ -A- Cg) d(x/d) 

Jo u 


From reference (a) 
1/d 


c - ft 

C Noc " i T 


j o ( '=^3 ) d(x/d > 


Hence * ^d =2 TTCf) _A- 


(4) 

(5) 


( 6 ) 


(7) 


( 3 ) 
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Thus* from eq # (8) and eq» (3) 



r/d and ^ are given in reference (c) as 


r/d = 


-A 2 - 


1/4) 2 - (x/d - f) 2 - (f 2 - 1/4) 


( LO) 


6 


tan 




f - x/d 


/d + - 1/4 


( 11 ) 


Hence f the running load over the tangent ogive may be evaluated by equation 
(9) utilizing the supplementary equations (4), (5)* (6), (10), (11) and 
appropriate charts and tables from references (b) and (d). 


The running load distribution over the cylindrical portion may be 
obtained via equation (2) and charts for the and G 2 functions from 
reference (b). Thus, the running load distribution for the ogive-cylinder 
is given by equations (9) and (2). 


CONCLUDING REMARKS 


In conclusion, the procedures and examples of application of the 
approximate second order shock expansion theory are presented herein to 
enable the user to predict the running load over a body of revolution 
at supersonic speeds# The process is quite simple and requires only a 
few hours of hand computations to establish the running load distribution 
and corresponding static stability characteristics# 


Edward E. Mayo 


Enclosures 
EEM: skd 
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SiMBOLS 


Ab 

C N*< 



d 

d c 

f 

Gi, G 2 , Gg, G^ 


body base area 

normal force coefficient curve slope, per radian, 
at * = 0 

pitching moment coefficient curve elope, p»r radlin, 
at ,jC ■ 0 

body base diameter 

cone diameter, d c (x/d) 

fineness ratio 
constants, reference (b) 


1 

M 

n 

P 

P t 

q 

r 

x 

x 

X 

V 

-A. 

a. 


body length 
Mach number 
running load, n (x/d) 
static pressure 
total pressure 

free- stream dynamic pressure 
body radius, r(x/d) 

longitudinal distance measured from nose 

longitudinal distance measured from nose-cylinder 
juncture 

angle-of-attack, radians 
ratio of specific heats 
surface slope, £ (x/d) 

loading, ^A- (x/d) 

Mach angle 

Prandtl-Meyer expansion angle 
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V 


SYMBOLS (Continued) 


lT~ 


cone half angle 
function defined by eq # 6 

ratio of cross-sectional area of streamtube to 
that at M = 1 


a 

c 

cp 

n 

s 

V 

tcv 


SUBSCRIPTS 

free-stream conditions 
afterbody 

quantities evaluated for cone tangent to body 

center of pressure 

nose 

quantities evaluated by generalized shock 
expansion method 

quantities evaluated at vertex of body 

quantities evaluated for cones tangent to body 
vertex 
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APPENDIX A 


RUNNING LOAD FOR A CONE- CYLINDER 
CONFIGURATION AT M = 6 


Consider the case of a 7.5° half angle cone configuration as shown 
in figure 1. From equation (1)* the running load value at the cone- 
cylinder juncture is 



/ r tan 



c 


For M = 6 t = 7.5°, from figure 4, reference (b) 

C N ^ = 1.875 
c 


Hence, 


= (3. 14) (. 1317) 0.875) = .774 

* qd 


From equation (2), the running load along the afterbody is given by 

7T 


qd j a 4 G 1 G 2 ex P[ " G 2 (fn + , 


where fn 


1 


= 3.8 


2 tan 2tan 7.5° 2(. 1317) 

and from figure 12, reference (b), (M = 6, ^ = 7.5°) 

G ? — 0.22, G -| = — *-52 — — — • 52 _ 

2 1 sin3> v .1305 " 3a 97 


Hence, 


n I Tr _ 

J = T 0.97) (.27) exp ^ - 0.22(3.8 + x/d) j 


= .686 exp ^ - .836 - 0.22 x/d j 
Thus, the afterbody running load values are 
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APPENDIX A 


x/ d 

x/d 

n 


x/d 

x/d 

n 

^ qd 

qd 

0 

3.8 

.297 


4 

7.8 

.123 

.5 

4.3 

. 266 


5 

8.8 

.099 

1.0 

4.8 

.239 


6 

9.8 

.080 

1.5 

5.3 

.214 


7 

10.8 

.064 

2 

5.8 

.191 


8 

11.8 

.051 

3 

6.8 

.154 



! 





The running load distribution is plotted in figure 1. The normal force 
curve slope, pitching moment, C m j and center-of-pressure location, 

x cp /d; are determined from the running load distribution as follows: 


C 


N 



(— — ) d(x/d) 
* qd 


(Al) 


C 


m 


_4 Ud 

rr ) 


(— n ) (x/d) d(x/d) 
qd 


( A2) 


x cp _ jc 

d ~ C N (A3) 


Integration of the load distribution in figure 1 for 1/d = ll e 8 (fa = 8) 
according to equations (Al), (A2) and (A3) yields 

^ = 3.28 per radian 


C.i, - 14.22 per radian 
x cp /d = 4.33 
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APPENDIX B 


RUNNING LOAD FOR AN OGIVE-CYLINDER 
CONFIGURATION AT M = 6 


Consider the case of a 3 to 1 fineness ratio nose ogive-cylinder 
configuration as shown in figure 2. The load distribution for the ogive 
is given by eq* (9) as 


" 2 J g 3 tan J v C Itr/ + G 4 tan ^ 

^ I 


+ Ul - G 3 ) C No< j ' G 4 j tan2 ^ 

L- / tcv -J '' 

t an ^ v 

I <• 

From eq. II \ = tan“l i f | = tan” (.343) = 18.93° 

£ f 2 - 1/4J 


From figures 3 and 4, reference (b), *^v 


Moo 

So, 


= .649; = 6.8 


= 1.774 


tcv 


Hence, M v = 6C.649) = 3.894 and from reference (d) 
P 

- .007635, V v = 64.302, ^ = 14.90, 


- 0 v 


= 9.7 


i/ 8 - & v + i/ v = 18.93 + 64.3 = 83.23 


M s a = 5 * 76 » P s a /P t = ‘O 0081 ^ . ^ = 9.998, 


. d - 44 - 7 


P. p s p t 

— - = ( — -)(— ) = 
p v s t 


.000813 

.007635 


=.1065 
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APPENDIX B 


P P ' P 

— =. | = (.1065) (6*8) = .724 

Poo ‘ 


i P v/ P ~/ 


Go is given by equation (4) as 


P_ sin 2 ^ 


G - 


3 P„ sin 2 


v _ e -Ta 


where y- a is given by equation (6) as 




/ (P s /P, ) M g z 
a a 


•o b v 


2(1 


P v 


a ) (M_ 2 - 1) 


( 1 . 4) ( .724) (5.76) 


9.7 _ 


. _ \ t / £- » / N J 1 f / V/ ^ / 

ra “ 2(1 - .724) (5.76) 2 - 1 44.7 


= 0.413 


Hence, 


G» = (.1065) (-^^>(.662) = .1024 
3 .342 

is given by equation (5) o*2> 

= 2(1 - e"^ a ) 


= 2(1 - . 662 ) = .676 

Thus, substitution into eq. (9) yields 

+ .676 tan ^ + 2.67 tan 2 & (Bl) 

by equations (10) and (11) as functions of the 


■ 2IT^.0623 
r 

where ^ and ^ are given 
fineness ratio and x/d« 
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APPENDIX B 


The solution of (Bl) yields 


x/d 

r/d 

tan ( 

n/ < qd 

0 

0 

.343 

| o 

.5 

.16 

.280 

i .462 

.7 

.21 

.257 

; .544 

1.0 

.28 

.221 

] .602 

1.2 

.32 

.198 

.605 

1.5 

.38 

. 164 

i .585 

2.0 

.45 

.109 

.474 

2.5 

.49 

.054 

i .327 

3.0 

.50 

0 

: .196 | 
J 1 


The above table gives the running load values over the ogive nose portion. 
The running load values along the afterbody are given by equation (2) as 


—I 

x q d l a 


lT 

TT 


G 1 G 2 


exp l -G 2 


Cfn + 


f) 

d J 


where fh . i= 3 and from figure 12, reference (b) (M = 6, ^ - 18.93°) 


Gj — .28; G, — 


.66 
sin £> , 


= 2.034 


Hence, 


= “£'(2.034) (. 28) exp £ (-.28) (3 + -x/d) 

= *447 exp / - .84 - .28 x/d , 

Thus* the afterbody running load values are 
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APPENDIX B 


x/d 

x/d 

1 

n/»<qd 

x/d 

x/d 

1 1 

n/u^gd 

X/d 

x/d 

1 

n^qd 

0 

3 

.193 ; 

5.0 

8 

.048 

10.0 

13 

.012 

.5 

3.5 

.168 j 

: 5.5 

8.5 

.041 

10.5 

13.5 

.010 

1.0 

4 

.146 : 

! 6,0 

9 

.036 

11.0 

14 

.0089 

1.5 

4.5 

.127 

i 6.5 

9.5 

.031 

11.5 

14.5 

.0077 

2.0 

5 

.110 

7.0 

10 

.028 

12.0 ! 

15 

.0068 

2.5 

5.5 

! .096 

; 7.5 

10.5 

.024 

12.5 i 

: 15.5 

.0059 

3.0 

6 

j .083 

; 8.o 

! li 

.021 ; 

13.0 

16 

.0049 

3.5 

6.5 

! .072 

8.5 

11.5 

.018 

14 

; 17 

j .0037 

1 4.0 

7 

.063 

9.0 

12 

.016 i 

15 

18 

1 .0028 

4.5 

7.5 

.055 

9.5 

12.5 

.014 : | 

16 

; 19 

: .0022 


x/d 

x/ d | 
1 

n/<*cqd 


I 

J 


17 

20 

.0017 

18 

21 

.0013 

19 ! 

22 

.0009 

20 

23 

.0007 

21 ! 

24 

.0005 

22 i 

25 

.0004 

23 | 

26 

.0003 

24 ' 

27 

.0002 

! 25 | 

1 

28 

.0002 
4 


The running load distribution is plotted in figure 2. Integration 
of the load distribution in figure 2 for 1/d = 13 (f a = 10) according 
to equations (Al), (A2) and (A3) yields 

=2.48 per radian 


Cm =7.22 per radian 


x /d = 2.91 
cp 


116 






118 






